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FLUTTER TESTS OF SOME SIMPLE MODELS AT A 
MACH NUMBER OF 7.2 IN HELIUM FLOW* 

By Homer G. Morgan and Robert W. Miller 


SUMMARY 


Results of hypersonic flutter tests on some simple models are pre- 
sented. The models had rectangular plan forms of panel aspect ratio 1.0, 
no sweepback, and bending-to -torsion frequency ratios of about 1 / 3 . Two 
airfoil sections were included in the tests; double wedges of 10 -, 
and 15 -percent thickness and flat plates with straight, parallel sides 
and beveled leading and trailing edges . The models were supported by a 
cantilevered shaft. 

The double -wedge wings were tested in helium at a Mach number of 7-2. 
An effect of airfoil thickness on flutter speed was found, thicker wings 
requiring more stiffness to avoid flutter. A few tests in air at a Mach 
number of 6.9 showed the same thickness effect and also indicated that 
tests in helium would predict conservative flutter boundaries in air. 

The data in air and helium seemed to be correlated by piston-theory cal- 
culations . Piston-theory calculations agreed well with experiment for 
the thinner models but began to deviate as the thickness parameter Mt 
approached and exceeded 1 . 0 . 


A few tests on flat-plate models with various elastic-axis locations 
were made. Piston-theory calculations would not satisfactorily predict 
the flutter of these models, probably because of their blunt leading 
edges. 

GuJtMrv 


INTRODUCTION 


Piston theory has been proposed by Lighthill (ref. l) as a simple 
method for calculating oscillatory aerodynamic forces at high Mach num- 
bers . Various other authors (for example, refs. 2 to 4) have developed 
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piston theory ahtL 'd^pliecf ft tt> ttertd'sttrtli^s Of* f lVftter characteristics. 
However, experimental data have not been available to determine the 
limits of applicability of this theory. Information which needs to be 
established includes the upper and lower Mach number bounds and the air- 
foil shapes and maximum airfoil thickness permitted by the theory. In 
addition, a predicted effect of airfoil thickness and shape on flutter 
speed has not been adequately demonstrated. 


In the present investigation, an attempt was made to isolate some 
of the thickness and shape effects predicted by piston theory by using 
models with double -wedge airfoils of varying thickness and some 
beveled-edge flat-plate models. The double-wedge models also permitted 
a partial evaluation of a restriction on piston theory given in refer- 
ence 1 - that w/a must be less than 1.0, where w is the downwash 
velocity and a is the free-stream speed of sound. For double-wedge 
models, this restriction reduces to Ut/a = Mt < 1.0 where U is the 
free-stream velocity, t is the thickness ratio and M is the Mach 
number. In order to assess these aerodynamic effects, the models chosen 
for the investigation were very simple to allow the structural portion 
of the flutter problem to be treated without difficulty. Also, the 
bending-to-torsion frequency ratio was kept low l/3) since 

trend studies (ref. 4) indicate that the largest effects of airfoil 
shape and thickness are in this range . 


The present investigation utilized helium as a testing medium to 
obtain flutter data at a Mach number of 7*2. Wind tunnels utilizing 
helium have previously been used extensively for studying hypersonic -flow 
phenomena. (See, for example, refs. 5 and 6.) These studies, however, 
have been conducted to examine steady flow problems such as total loads 
or local pressure distributions. Meanwhile, very little flutter data 
exist at hypersonic speeds. Most of the available information is given 
in references 7 and 8 and was obtained at a Mach number of 6.9 in air. 

Helium has two primary advantages as a hypersonic test medium when 
compared with air: (l) Higher Mach numbers can be obtained without 

requiring the addition of heat to a system since helium has a much lower 
liquefying temperature than air. (2) Higher dynamic pressures can be 
obtained at the same Mach number and stagnation pressure because helium 
has a larger ratio of specific heats than air. Thus, use of helium 
results in a less complicated and less expensive wind-tunnel facility to 
produce the same test conditions. The higher dynamic pressures available 
in helium as compared with air at the same storage pressure and the feasi- 
bility of using atmospheric stagnation temperatures with helium are 
especially attractive for flutter testing since model design and con- 
struction are greatly simplified. 
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A few test runs were made in air at a Mach number of 6.9 in an 
attempt to correlate data on flutter models tested in two mediums . 


SYMBOLS 

a free-stream speed of sound, ft/sec 

b semichord, in. 

Cp pressure coefficient 

c chord, in. 

d shaft width, in. 

f frequency, cps 

I a mass moment of inertia about the elastic axis per unit 

length, slug-ft^/ft 

M Mach number 

m mass per unit length, slugs /ft 

P Q stagnation pressure, Ib/sq ft 

q dynamic pressure, lb/sq ft 

boi, . — 

R — yj[i, a stiffness-altitude parameter 

static unbalance about the elastic axis per unit length, 
positive for center of gravity back of elastic axis, 
slug-ft/ft 

t shaft thickness, in. 

U free-stream velocity, ft/sec 

w downwash velocity, ft/sec 

x Q elastic-axis position in fraction of the chord, measured from 

the leading edge 


• • 


tt in i ..14,#%;* w t 
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H mass of moS.e*L dfvide*(f by*mas*s *of fluicf ‘contained in a cylinder 

determined by the model, m/itpb 2 

p density of test medium, slugs/cu ft 

t maximum thickness of model in fraction of the chord 

a) circular frequency, radians/sec 

Subscripts : 

c calculated quantity 


e 

f 

h 

a 


experimentally determined quantity 

quantity at flutter 

bending mode 

torsion mode 

free -stream conditions 


APPARATUS 


Langley Hypersonic Aeroelasticity Tunnel 


Most of the tests were performed in the Langley hypersonic 
aeroelasticity tunnel which uses helium as a test medium. A sketch of 
this blowdown tunnel is shown in figure 1. The nozzle used in the tests 
has circular cross sections with a 1 . 4l-inch-diame ter minimum section 
and an 8-inch-diameter test section. Helium is supplied to the stagna- 
tion chamber at pressures up to 1,200 pounds per square inch frcan which 
dynamic pressures of over 5,000 pounds per square foot are obtainable. 
The downstream end of the tunnel is connected to a vacuum chamber which 
can be operated at pressures as low as l/2 inch of mercury absolute. 
Stagnation temperature is essentially constant during a run and corre- 
sponds to atmospheric temperature. With the available high-pressure 
helium supply, test runs are of approximately 50 -second duration. 

Models are wall-mounted on a removable plate which is contoured to 
maintain a circular test section. The opposite wall is made of trans- 
parent plastic; thus, the model could be observed during the test runs. 
Mach number surveys have been made across the diameter of the tunnel at 
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the center of the test section and 4 inches upstream and li inches down- 

o 


stream from the center of the test section. The results are shown in 
figure 2. The Mach number at the center of the test section, where the 
models are located, is about 7.2 and remains practically constant over 
the span of the model. 


Langley 11 -Inch Hypersonic Tunnel 

Three test runs were made in the Langley 11-inch hypersonic tunnel 
which produces a Mach number of 6.9 when air is used as a test medium. 
The characteristics of this tunnel are given in reference 9- 


MODELS 


Configuration 

All models tested had panel aspect ratios of 1.0, no taper, zero 
sweep, and 4 -inch span. The models were supported by a s haf t, which was 
2 inches long. The shaft was cantilevered at the tunnel wall; thus, the 
model was held outside the boundary layer and in a region of uniform 
flow. Strain gages for determining model response were attached to the 
shaft. The shaft was covered by a fairing which was rigidly attached to 
the tunnel wall. Provision was made for a clamping device which extended 
from the tunnel wall along the leading edge of the fairing and restrained 
the model during the starting transient. However, this device was found 
to be unnecessary and its use was discontinued after the first few tests. 
A photograph of a typical model, ready for installation in the test sec- 
tion, is shown in figure 3- 


Construction 

Details of the model construction are shown in figures 4(a) and 4(b) . 
Two types of models were tested: models 1 to 4 were flat plates with 

beveled leading and trailing edges (fig. 4(a)) whereas models 5 to 17 
were double-wedge airfoils (fig. 4(b)). 

The flat-plate models were solid aluminum alloy with an integral 
shaft. The center line of the shaft was varied along the chord of the 
wing to change the elastic-axis position. The leading and trailing edges 
were beveled at a 45° angle. 

The double -wedge models had a steel core which was integral with the 
shaft. Balsa, which formed the airfoil contour, was glued to the core 
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with its grain oriented perpendicular to the core to minimize the stiff- 
ness contribution of the wood. Part of the core was drilled out, before 
covering with balsa, to reduce the model weight. Models of 5~> 10-, and 
15-percent thickness were tested. The elastic axis of this group of 
models was at 45 percent chord. 


Physical Parameters 

The mass parameters of the models are listed in table I along with 
pertinent dimensions. These data were measured prior to installing the 
models in the tunnel. The mass of the shaft is not included in the 
data shown. 


All models were vibrated with an acoustic shaker before testing 
and the vibration nodes and frequencies were determined. Typical nodal 
patterns for both types of models sire shown in figure 5- lu all cases 
examined, the third and fourth frequencies, corresponding to second 
bending and second torsion, were well above the first and second fre- 
quencies. ^As an example, model 10 had frequencies of fj, = 20.5 C P S , 

f =64.5 cps, fj^ = 188 cps, and = 402 cps.^ Since the flutter 


frequency always fell between first bending and first torsion, the third 
and fourth frequencies were not determined for all cases. The natural 
frequencies are listed in table I. 


TEST PROCEDURE 


Models were mounted in the test section at zero angle of attack. 
Then, after installation in the tunnel and just prior to the test run, 
the first two natural frequencies of the model were checked. 

For the flutter tests, flow was established in the tunnel at a low 
dynamic pressure. Then, with the Mach number remaining constant, dynamic 
pressure was increased until the model fluttered or maximum operating 
conditions were reached. Stagnation temperature and pressure and signals 
from the strain gages mounted on the model shafts were recorded on an 
oscillograph throughout the test. Conditions at flutter were established 
by correlating the stagnation quantities recorded on the oscillograph 
record with the model response given by the strain -gage signals. 
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Flutter calculations were made on the models by using two 
uncoupled modes. The modes used were chosen after examination of pho- 
tographs of the modes for similar models. These photographs indicated 
that it would be reasonable to assume that the torsion mode consisted 
of unit twist along the entire span and the bending mode was the 
straight line given by 


F h = l/Mi + 3n) 

(F^ is the displacement at the elastic axis and T) is a spanwise 
variable which is zero at the root and unity at the tip.) The calcu- 
lated flutter results were found to be relatively insensitive to the 
slope of this line . The frequencies associated with these modes were 
those which were measured and axe listed in table I. 

The distributed mass properties used in the calculations were 
those listed in table I. The mass of the shaft is not included in these 
calculations . 

The aerodynamic forces were obtained from second-order piston 
theory. The coefficients axe presented in reference 2. Airfoil-shape 
terms were included for the flat-plate models as well as for the double - 
wedge models. 


RESULTS AND DISCUSSION 


The results of the tests axe listed in table II. The test medium, 
frequency ratios, and test-section conditions are shown at flutter. For 
tests where no flutter occurred, the conditions which axe given corre- 
spond to the maximum dynamic pressure reached during the test. A 


stiffness -altitude parameter 



is also shown at flutter or at 


the maximum dynamic pressure. This parameter has been shown by piston 
theory calculations to be useful for correlating flutter data at high 
Mach numbers for the range of density ratios covered by the investigation. 


The stiffness -altitude parameter depends only upon the physical 
properties of the wing - in particular, the torsional stiffness - and 
upon the atmosphere in which it operates . Its value increases as either 
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the altitude or the stiffness increase. The stable region, when this 
parameter is the ordinate, will be above the flutter boundary. 


Double -Wedge Models 

Experiment . - The experimental results from table II for the double - 
wedge models are plotted in figure 6 as functions of the thickness 
parameter Mr. As mentioned earlier, the parameter Mt has certain 
significance concerning the limits of piston-theory applicability. The 
ratio of flutter frequency to torsion frequency is shown in figure 6 (a). 
This frequency ratio exhibits considerable scatter at any particular 
value of Mr but is relatively constant with thickness (or Mr). The 
flutter frequency varies from 4 5 percent to 65 percent of the torsional 
frequency, whereas the bending frequency is about one-third the tor- 
sional frequency. 

Examination of the stiffness-altitude parameter for the double 
wedges (fig. 6 (b)) reveals a definite thickness effect. As Mr 
increases, the value of the stiffness -altitude parameter increases; 
thus, flutter occurred at a higher altitude or with greater stiffness. 
Although considerable scatter is evident from the data, the general 
trend is for the thickness effect to be nonlinear with thickness and 
to decrease as Mt increases. 

Also shown in figure 6 are data obtained both in air (at M = 6.9) 
and in helium (at M = 7 . 2 ). Only three runs were made in air; the 10 - 
and 15-percent-thick models fluttered and the 5 - P er cent- thick model 
would not flutter at the available dynamic pressure. The flutter fre- 
quency, obtained from the models which did flutter, is seen from fig- 
ure 6(a) to correlate very well with flutter frequencies determined in 
helium. From figure 6(b), the stiffness -altitude parameter at flutter 
is seen to be 20 to 25 percent lower in air than in helium on the two 
models which fluttered. This result would indicate that studies made 
in helium at this Mach number would predict flutter boundaries which 
would be conservative in air. However, as will be pointed out later, 
these differences may be accounted for by using piston-theory 
calculations . 

Calculations . - The thickness effects predicted by piston theory, 
which is used to determine aerodynamic forces in these calculations, 
change drastically with frequency ratio. The models in the test pro- 
gram were designed to have bending-to-torsion frequency ratios of 
approximately one-third, in the range where the predicted effects are 
large. A comparison between calculations and experiment on models 5 , 

6, and 7 is shown in figure 7 to illustrate this variation with 
frequency ratio. Figure 7 shows that the calculations agree very well 
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with experiment on the 5- and- 10-percent-thick models which were tested 
at M = 7.2 in helium. Calculations on the 15-percent-thick wing pre- 
dict a greater effect of thickness than is measured. As the frequency 
ratio approaches 1.0, the calculations show flutter at much higher 
altitudes, although the magnitude of the thickness effect is much 
reduced. 

A summary of flutter calculations on all the double-wedge wings 
plotted against the thickness parameter Mr is shown in figure 8. The 
ratio of experimentally determined flutter frequency to calculated 
flutter frequency is given in figure 8(a). In all cases, calculated 
flutter frequencies are too high. However, air and helium data all fall 
within the same scatter band. 

The ratio of experimental stiffness -altitude parameter to calcu- 
lated stiffness-altitude parameter is plotted in figure 8(b). The trend 
is for the calculations to give poorer agreement with experiment as the 
thickness parameter Mr becomes larger; thus piston theory is giving a 
poorer description of the aerodynamic forces. At fir = O.56, calcula- 
tions are within 8 percent of experiment; at Mt « 0.70, calculations 
are all within l6 percent of experiment; and at Mr « 1 . 05 , calculations 
are all within 25 percent of experiment. A check calculation on some 
models using third-order piston theory actually showed poorer agreement 
with experiment than the calculations with second-order piston theory. 

As shown in figure 8(b), flutter points obtained in air are calcu- 
lated with about the same degree of accuracy as those obtained in helium. 
Calculations in air also show that the increase in stiffness -altitude 
parameter when the Mach number is changed from 6.9 to 7*2 is only about 
4 percent. This difference is not as great as the experimental differ- 
ence between air and helium shown in figure 6(b) . However, the calcula- 
tions in figure 8(b), which contain the variation in 7, the ratio of 
specific heats, do tend to correlate the data and suggest that piston- 
theory calculations may be used to extrapolate flutter data obtained 
in helium to air. The flutter calculations are affected by a change in 
the ratio of specific heats since this factor 7 appears in the second 
term of the second-order piston-theory pressure coefficient 



The thickness effect arises from the second term of this relation and, 
since 7 is 1.67 for helium compared with 1.4 for air, flutter calcu- 
lations in helium will show a larger thickness effect than flutter cal- 
culations in air. 
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As mentioned previously, a restriction on the use of piston theory, 
as given in reference 1, is that w/a must be less than 1.0. If the 
downwash due to unsteady motion is neglected, this requirement for 
double-vedge airfoils becomes Mt < 1.0. (This condition is reasonable 
since flutter calculations are made for small perturbations at the 
neutral stability point.) This limitation on piston theory seems to be 
borne out at M = 7.2 by the comparison shown in figure 8. 


Flat-Plate Models 

Experiment . - The experimental results for the flat-plate models are 
shown in figure 9 where the stiffness-altitude parameter is plotted as 
a function of elastic -axis location in fraction of the chord. Figure 9( a ) 
shows that the flutter frequencies are all less than 40 percent of the 
torsion frequency. Thus, flutter frequencies were about the same as the 
bending frequencies, which were approximately one -third of the torsion 
frequencies. By observation, the flutter mode was determined to be 
almost pure bending with only slight torsion components. 

The stiffness-altitude parameter at flutter is plotted in fig- 
ure 9(b) . The no-flutter point was found on a model (number l) whose 
frequencies were about 65 percent higher than those of the model with 
the same elastic-axis location which did flutter at a lower dynamic 
pressure (model number 2). The difference in the stiffness-altitude 
parameters possibly is explained by noting that the center of gravity 
of the model which fluttered is slightly farther back of the elastic 
axis whereas its radius of gyration is slightly less than the model 
which did not flutter. 

Calculations . - A comparison of experiment and piston-theory calcu- 
lations plotted against elastic-axis location is shown in figure 10. 

The ratio of experimental to calculated flutter frequencies is presented 
in figure 10(a) . The experimental flutter frequencies are approximately 
one-half the calculated flutter frequencies. In figure 10(b), the ratio 
of experimental to calculated stiffness -altitude parameter is shown. 

Poor agreement of experiment and calculations is apparent for forward- 
axis locations, the calculated value being only 52 percent of the experi- 
mental value at x Q = O.388. The agreement improves as the axis moves 

rearward until the calculated value is about 95 percent of the experi- 
mental value at x Q = 0.513. 

These flat-plate models have airfoil sections which are straight- 
sided with a bevel at both the leading and trailing edges, as shown in 

figure 4. The bevel amounts to only 1^ percent of the chord at each 
end of the model. This bevel appears to be such a small fraction of 
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the chord that the airfoil would be expected to behave like a zero- 
thickness plate. (For flat sides, the slope of the airfoil contour is 
zero and second-order piston theory will not give a thickness effect.) 
However, all calculations made on these models by assuming zero thick- 
ness gave infinite flutter speeds. The calculations presented in 
figure 10 use the piston-theory thickness term, which includes the 
effect of the bevels at the edges, despite the fact that the requirement 
for w/a < 1.0 is violated in the region of the bevel. Using this term 
apparently overestimates the shape effect, at least for forward-axis 
locations . 

A possibility for explaining the poor agreement of experiment and 
calculations is the inability of the theory to account for a blunt 
leading edge. Studies of blunt-nosed slender bodies (for example, 
ref. 10) have revealed a pressure carryover along the forward portion 
of the body due to the blunt leading edge and its detached shock wave. 
Since flutter of an aerodynamic surface is known to be very sensitive 
to the relative position of the center of pressure and elastic axis, 
such a pressure carryover on an oscillating airfoil could well have 
large effects on the flutter characteristics. 


CONCLUDING REMARKS 


Tests at a Mach number of 7*2 in helium on simple models having 
double-wedge airfoils of various thicknesses indicate a large effect of 
thickness on flutter speed. As airfoil thickness increases, the models 
require more stiffness at a given altitude to avoid flutter. The 
thickness effect on flutter speed is nonlinear, its magnitude decreasing 
with increasing airfoil thickness. 

A limited number of tests on similar models at a Mach number of 6.9 
in air indicate the same trends with thickness as were observed in 
helium tests. Comparison of tests in air and helium indicates that 
flutter boundaries obtained in helium would be conservative in air, 
possibly because of the differences in the ratios of specific heat. 

Piston-theory calculations on the double -wedge models agree fairly 
well with experiment for the thinner airfoils . As the thickness 
increases, such that the thickness parameter Mr approaches and exceeds 
one, agreement between the calculations and experiment becomes poorer. 
Tests in air and helium were calculated with about the same accuracy and 
are correlated by the piston-theory calculations; thus tests in helium 
can be extrapolated to air conditions with some confidence. 

The flutter of some flat-plate models having beveled edges could not 
be calculated with reasonable accuracy by using the piston theory. The 



failure of piston theory to calculate flutter on these models can prob- 
ably be attributed to the effect of the very blunt leading edges of the 
model. 


Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Field, Va., January 9 > 1959 - 
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Figure 3.- Photograph of a typical model ready for installation in the L- 58 - 3035 

Langley hypersonic aeroelasticity tunnel. 
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Second mode 

Third mode 

Fourth mode 



(a) Flat -plate models. 



(b) Double-wedge models. 

Figure 5.- Typical no dal patterns for the first four vibration modes of 

the test specimens. 
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(a) Ratio of flutter frequency to torsional frequency 
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(b) Stiffness -altitude parameter. 

Figure 6.- Experimental flutter results on double -wedge models plotted 
against the thickness parameter Mt. 0.31 = 0.34; Xq = 0.4 
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Figure 7 .- St iff ness -altitude parameter at flutter plotted against fre- 
quency ratio for double -wedge models tested in helium at M = 7-2. 
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(a) Ratio of experimental flutter frequency to calculated 

flutter frequency. 
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(b) Ratio of experimental stiffness-altitude parameter to calculated 

stiffness-altitude parameter. 

Figure 8.- Comparison of experimental and calculated flutter results for 
double-wedge models at various values of the thickness parameter Mr . 

0.31 ^ x o = R = Vp"* 
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(a) Ratio of flutter frequency to torsional frequency. 
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(b) Stiffness -altitude parameter. 

Figure 9-- Experimental flutter results for flat-plate models plotted 
against elastic-axis location, o^/o ^ = 0.33 and. 0.3^-j M = "J .2 in 

helium; R = 


64t- 



(a) Ratio of experimental flutter frequency to calculated flutter 

frequency . 
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(b) Ratio of experimental stiffness-altitude parameter to calculated 

stiffness-altitude parameter. 

Figure 10.- Comparison of experimental and calculated flutter results 
for flat-plate wings at various elastic-axis locations. = 0. 

and O.jk; R = J\T. 
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